Microwave Background Signals from Tangled Magnetic Fields 
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An inhomogeneous cosmological magnetic field will create Alfven-wave modes that induce a small 
rotational velocity perturbation on the last scattering surface of the microwave background radiation. 
The Alfven-wave mode survives Silk damping on much smaller scales than the compressional modes. 
This, in combination with its rotational nature, ensures that there will be no sharp cut-off in 
anisotropy on arc-minute scales. We estimate that a magnetic field which redshifts to a present 
value of 3 x 10 -9 Gauss produces temperature anisotropies at the 10fj,K level at and below 10 arc- 
min scales. A tangled magnetic field, which is large enough to influence the formation of large scale 
structure is therefore potentially detectable by future observations. 



PACS Numbers : 98.80.Cq, 98.70.Vc, 98.80. Hw 

The origin of ordered, large-scale cosmic magnetic 
fields remains a challenging problem. It is widely as- 
sumed that magnetic fields in astronomical objects, like 
those in galaxies, grew by the amplification of small seed 
magnetic fields by turbulent dynamo action pj. How- 
ever, the efficiency of the turbulent galactic dynamo is 
still being debated [Q. Alternatively, the galactic field 
could be a remnant of a cosmological field of primordial 
origin ||, although, as yet, there is no entirely compelling 
mechanism for producing the required field. It could be 
present in the initial conditions, be produced quantum 
gravitationally or at a phase transition, or be generated 
in some way at the end of a period of inflation Q] . A pri- 
mordial field that expanded to contribute a present field 
strength of order 10~ 9 Gauss, tangled on galactic scales, 
could also impact significantly on galaxy formation . 
It is of considerable interest, therefore, to find different 
ways of limiting or detecting such primordial fields 
We argue that observations of anisotropies in the cos- 
mic microwave background radiation (CMBR), provide 
a potentially powerful constraint on such fields. Indeed, 
earlier work has found that the isotropy of the CMBR 
already places a limit of 6.8 x 10~ 9 (il h 2 ) 3 Gauss on the 
present strength of any uniform (spatially homogeneous) 
component of the magnetic field |!|. (fig is the present 
density parameter, and h the Hubble constant in units 
of lOOkms Mpc .) Here, we obtain comparable con- 
straints on tangled (inhomogeneous) magnetic fields as 
well and highlight the distinctive fluctuation signature 
that they are expected to leave in the small-scale struc- 
ture of the CMBR. 

We build on the results of an earlier paper || , where we 
considered the evolution and damping of cosmic magnetic 
fields in the early universe. We showed that inhomoge- 
neous magnetic fields induce rotational Alfven modes as 
well as compressional density and velocity perturbations 
(see also ref. JtJ ) . The compressional fluid perturbations 
are of small amplitude (~ V\ << 1) on all scales, be- 
cause of the restoring force contributed by the pressure 
of the radiation-baryon fluid. Here, Va is the Alfven ve- 
locity in the radiation era of the universe (see eqn. (|^) 



below) and we set c = 1. Compressional modes are also 
strongly damped on scales smaller than the comoving 
photon diffusion length, Ls, by Silk damping (To). How- 
ever, the rotational component of the fluid velocity, in- 
duced by tangled magnetic fields, can be much larger (a 
significant fraction of Va)- Two features are important. 
First, the restoring force is much weaker (since radiation- 
baryon pressure does not affect this mode), and further- 
more, these velocity perturbations survive damping on 
comoving scales larger than La ~ VaLs @- It is these 
rotational velocity perturbations, induced by the tangled 
field, which can produce significant small angular scale 
anisotropies in the CMBR, through the Doppler effect. 

We calculate this effect in two steps. First we calcu- 
late the anisotropy at small angular scales produced by a 
general rotational component of the velocity, due to the 
Doppler effect. Then we estimate the magnitude of this 
rotational component for a given spectrum of magnetic 
inhomogeneities. 

The equation governing the evolution of temperature 
perturbation AT(x, 7, tq)/T$, can be derived from the 
moments of the Boltzmann equation for photons [ pd| . 
Here, x is the comoving position, r is conformal time, To 
its present value, and 7 is the direction of photon prop- 
agation. It is conventional to expand AT/Tq in spheri- 
cal harmonics with AT/Tq — J2i m a imYi m {l) and define 
C; =< |a/ m | 2 >. For rotational velocity perturbations a 
detailed calculation of Ci has been done in Ref. and 
p3[ in the context of computing the Vishniac effect, and 
we can adopt their final result. We have (from Eqs. (5) 
and (16) of pi), 
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Here, V(k, r) is the Fourier transform of the rotational 
component of the fluid velocity, with k the co-moving 
wave number, and ji(z) is the spherical Bessel func- 
tion of order I. The 'visibility function', <?(to,t), de- 
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termines the probability that a photon reaches us at 
the epoch tq if it was last scattered at the epoch r. 
It is given by g(r,r') = reexp[— / , k(t")cIt"], where 
k(t) = rt e (r)<7Ta(T), n e is the electron number density, 
<7t is the Thomson cross section, and a(r) is the cosmo- 
logical scale factor normalised to unity at the present. 
We also adopt a flat universe throughout. 

For standard recombination physics, g is peaked about 
a small range of conformal times around the time of re- 
combination. For the present purpose of making approxi- 
mate analytical estimates, it suffices to approximate it as 
a Gaussian, and take g(To,r') = (2tt<7 2 )~ 1 / 2 exp[— (r' — 
t») 2 / (2a 2 )], where t* is the conformal epoch of "last scat- 
tering" and a measures the width of the last scattering 
surface (LSS). Using the expressions given in Ref. juj, 
we estimate r* ~ 178.2h~ 1 Mpc and a = ll,lh~ x Mpc. 
Note that due to the visibility function, the dominant 
contributions to the integral over r in Eq. (|l|) come from 
conformal times in a range a around the epoch r = r*. 
The presence of ji(k(ro — r)) picks out (fc,r) values in 
the integrand which have fc(r — r) ~ /. 

We can obtain analytic estimates of C\ in two limits. 
First, for wavelengths such that ka << 1, fc(r — r), 
and hence ji(k(jQ — r)), vary negligibly over the range 
of t where g is significant. So they can be evaluated at 
r = r* and taken out of the integral over r in Eq. (Q). 
Furthermore, in general, V(k,r) does not vary rapidly 
with conformal time within a, nor does it vary rapidly 
with k in the k range around k ~ l/R* where ji(kR*) 
contributes dominantly (we define = r — r*). Thus, V 
can also be evaluated at r = r* and k — l/R* and pulled 
out of the integrals. The remaining integral of g over r 
gives unity, while that over jf can be done analytically, 
giving 
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where, A^(fc,r*) = fc 3 < |F(fc,T»)| 2 > /(2tt 2 ) is the 
power per unit logarithmic interval of fc, residing in the 
rotational velocity perturbation V. 

In the other limit, ka >> 1, for wavelengths much 
smaller than the thickness of the LSS, one can treat g as 
a slowly-varying function compared to ji in the integral 
over r in (|l|) . There is a cancellation due to superposition 
of oscillating contributions of ji over the thickness of the 
LSS. An approximate evaluation of the angular power 
spectrum in this case gives 
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Note that in this small- wavelength case, (ka >> 1), the 
angular power-spectrum is suppressed by a l/ka factor 
due to the finite thickness of the LSS. We can now eval- 
uate the CMBR anisotropics produced by tangled mag- 
netic fields. 

We will assume that the tangled magnetic field is ini- 
tially a Gaussian random field. When the Hubble ra- 
dius grows to encompass tangles on any particular scale, 



the Lorentz force associated with the magnetic inhomo- 
geneity induces a fluid velocity. One would like to es- 
timate the rotational component of this velocity, for a 
general spectrum of magnetic inhomogeneities. Before 
doing this however, we first consider the simpler example 
of the special non-linear cosmological Alfven wave solu- 
tion discussed in jgj. This will allow us to estimate in a 
simpler context the magnitude of the CMBR anisotropy 
induced by a tangled magnetic field. 

We showed in ||, that the equations of relativistic, 
viscous, magnetohydrodynamics are conformally invari- 
ant during the epochs when the radiation energy density 
dominates over the baryon density. We exploited this to 
transform the MHD equations from the Friedmann uni- 
verse to flat spacetime, and studied a special nonlinear 
Alfven mode solution. In this solution, the magnetic field 
B is split into a uniform field, say along the z direction, 
and a tangled component of arbitrary strength, perpen- 
dicular to z. All variables depend only on z and r. In 
particular, B = [B z + b (z, r)n]/a 2 , where n.z = 0. 
The tangled field, bo(z,r), can be arbitrarily large com- 
pared to B and satisfies a linear damped wave equation. 
Expanding the magnetic inhomogeneity in plane wave 
modes, as bo(z,r) — f(k,r)e lkz , this is, 0, 
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Here, derivatives are with respect to conformal time. The 
associated velocity perturbation is given by v(z, r) = 
v o(k, T)e lfez n; from the induction equation, vo(k,r) = 
—if/(kB ). InEq.(||), rj = (4/15)p 7 / 7 is the shear viscos- 
ity coefficient associated with the damping due to pho- 
tons, whose energy density is p 7 and mean-free-path is 
Z 7 . We have also defined the Alfven velocity Va by 
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where p and p are the energy density and pressure of 
the relativistic radiation-baryon fluid; po and po their 
redshifted present-day values and B the present field 
strength assuming the field is simply frozen into the uni- 
formly expanding universe. To obtain a numerical esti- 
mate, we have expressed -Bo hi units of 10~ 9 Gauss and 
taken the dominant contribution to the energy density 
p = p~ n the photon energy density, as is appropriate in 
the radiation era after the epoch of e + e~ annihilation. 

We can solve Eq.(|]) analytically in two limits. For 
scales (or fc -1 ), larger than the comoving Silk scale 
at recombination, which is L$ ~ (Z 7 (i*)t*) 1 / 2 /a(t*), 
and assuming the fluid starts from rest, we can ne- 
glect the frictional term compared to the driving term 
in eq. (^). (Here t is the usual cosmic time). 
Since /(rj) = at some initial time, say Tj, we 
have / = /o(fc) cos[fcV^(r — t;)] and vo(k, t) = 
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(fo(k)/iB Q )VA sm[fcV^(r — r,*)]. For a flat universe, com- 
posed of radiation and matter, we can integrate the 
Friedmann equation to obtain r = 6000/i -1 [y / aTTE^ — 
Mpc. The phase of oscillation is then x — 
W a {t - n) « 10- 2 B_ 9 (fc/(0.2/iAfpc^ 1 ))(r/r„), where 
we assume r* >> Tj, and adopt a LSS redshift, z* ~ 
1100. Since X « 1, Wb(*,r*) « (/ (fe)/iB )feVf n - 
0.86 x 10- 5 (/ (fc)/i-B )(fc/0.05/iMpc' 1 )( J B_ 9 /3) 2 . (One 
can also check from this solution that the neglect of the 
/ term is a good approximation for kLs < 1). The asso- 
ciated CMBR anisotropy will depend on the orientation 
of the tangled and large-scale fields with respect to the 
line of sight. A rough estimate of the maximum signal 
for kLs << 1 is 

AT - Toluol ~ 23.3 f iK{k/0.05hMpc- 1 ){B- 9 /3) 2 ; (6) 

where To is the present mean temperature of the smooth 
CMBR. We have assumed that the tangled and uni- 
form field components are of the same order. (Also, on 
these scales ka < 0.5, and there is negligible damping of 
anisotropics due to the finite LSS thickness). 

In the opposite limit kLs » 1> damping due 
to photon viscosity dominates, and one can as- 
sume the fluid reaches a terminal velocity, where 
/ » 0, so v « (if/B )5V!/(kL (r)) ~ 0.76 x 
W~ 5 (if/B )(B_ g /3) 2 (k/0.25h)- 1 f b h- 1 . Here, the 
quantity L 1 = i 7 (r*)/a(r*) ~ 3AMpcf^~ , is the comov- 
ing photon mean-free-path at last scattering. We have 
also defined /& = (flb/0- 0125ft- -2 ), Of, being the present 
baryonic density parameter (see Eq. (4.12) of Ref. [p[) . 
The oscillator is overdamped on the relevant scales, [pi . 
Thus, the solution with /(n) = will have / ~ / . Fur- 
ther, since L5 and a are comparable, we have ka > 1, 
and so there will be a {ka) -1 / 2 suppression in the ob- 
served anisotropy due to the finite thickness of the LSS. 
Therefore, in this case we predict that, for kLs >> L 
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Again, we have assumed comparable tangled and homo- 
geneous fields. The geometrical factor Co < 1 incor- 
porates the unknown orientation of the large-scale field 
and the small-scale tangles, together with the precise 
form of the LSS thickness suppression effects. We see 
that in both limits, (||)-(@) small angular scale CMBR 
anisotropies ~ lOpK could arise for magnetic fields Bq 
- 3 x 10 -9 Gauss. 

Let us now turn to a more general calculation of the 
magnetically induced rotational- velocity perturbation, in 
the case when the fields are initially stochastic with Gaus- 
sian probability distribution. On galactic scales and 
above, the induced velocity is generally so small that it 
does not lead to any appreciable distortion of the ini- 
tial field (see above and ||). So, to a very good ap- 
proximation, one can assume that the evolution of the 



magnetic field is simply a dilution due to the uniform 
Hubble expansion, with B(x, t) = b (x)/a 2 . (Here we 
assume bo does not include any uniform component.) 
The Lorentz force associated with the tangled field is 
then given by Fl — (V x b ) x b /(47ra 5 ). The Lorentz 
force still pushes on the fluid, causing rotational velocity 
perturbations. We also focus on the perturbations with 
co-moving length scales larger than the photon mean- 
free-path at decoupling, and describe the viscous effect 
due to photons, as before, in the diffusion approximation. 
The Fourier transform of the linearised Euler equation for 
the rotational perturbations in the baryon-photon fluid 
is given by 
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Here we have defined the Fourier transforms of the mag- 
netic field, b (x) = J2k b(k) exp(ik.x) and the velocity 
field, v(x, t) — V(k, t) exp(ik.x), where V.k = for 
the transverse (vortical) part of the velocity. Also pb 
is the baryon density. The fact that the Lorentz force 
is non- linear in bo(x), leads to the mode-coupling term 
F(k) = E p [b(k + p).b*(p)]p - [k.b*(p)]b(k + p), and 
the projection tensor, -Py (k) = [Sij — kikj/k 2 ] projects F 
onto its transverse components (perpendicular to k ). 

As before, we can solve Eq.(0) analytically, on scales 
much larger or much smaller than the Silk scale. For 
kL s < 1, and when the fluid starts from rest (Vi(rj) = 0), 
the damping due to the photon viscosity can be ne- 
glected compared to the driving due to the Lorentz 
force. Integrating the resulting first-order equation gives 
Vi = Gi{r — Ti)/(1 + S*), where we have defined d — 
P^Fj/lAiripo + po)] and S* = (3p 6 /4p 7 )(n) ~ 0.4/ 6 . 
The (1 + jS*) factor results from the reduction in the in- 
duced velocity due to baryon inertia, an effect neglected 
in the non-linear Alfven wave solution above. One can 
also check from this solution that the neglect of viscous 
damping is valid for kLs < 1- In the other limit, with 
kL s >> 1, we can use the terminal- velocity approxima- 
tion, neglecting the inertial terms in the Euler equation, 
and simply balance the Lorentz force by friction. This 
leads to Vi = (G i /fc)(fcT 7 /5)- 1 . 

To calculate the induced CMBR anisotropy, we need 
to define the spectrum of the tangled magnetic field, 
M(k). We define, < 6 i (k)6 i (q) >= 5 k , q P^(k)M(fe), 
where <5k, q is the Kronecker delta which is non-zero 
only for k = q. With this definition we also have 
< b 2 >= 2f(dk/k)A 2 (k), where A 2 (fc) = k 3 M (k) / {2n 2 ) 
is the power per logarithmic interval in k space resid- 
ing in magnetic tangles, and we have changed the sum- 
mation over k space to an integration. The ensemble 
average < \V\ 2 > = < V t V* >, and hence the C;s, 
can be computed in terms of the magnetic spectrum 
M(k). It is convenient to define a dimensionless spec- 
trum, say h(k) — A 2 (/c)/(£? 2 /2), where Bq is a fiducial 
constant magnetic field. We will also define the Alfven 
velocity, Va, for this fiducial field, using Eq.(||). Fi- 
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nally, as a measure of the CMBR anisotropy induced 
by the tangled magnetic field, we define the quantity 
AT B (l) = {l(l+ l)Q/27r}V2T . 

Since scales with kL s < 1 also generally satisfy the 
criterion ka < 1, the resulting CMBR anisotropy on 
these scales can be estimated using Eq.(||). A lengthy but 
straightforward calculation gives, for scales with kL s < 1 
and kcr < 1, 



AT B (0 
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0.4, and we have again 



Here, I = as before, S"* 
assumed t* >> T{. 

For scales with khg > 1 and fcer > 1, we can use Eq.(||), 
and Vi = (Gi/k)(kL 7 /5)^ 1 . In this case a similar calcu- 
lation to the one performed above gives 



7T 1 / 4 

AT B (l)=T ^=I(k) 



B. 



w 2 



kL^r^ka) 1 / 2 



5.6/iif 



/ 



1500 



-3/2 



KoH^) (io) 



where /150 = (ft-/0.5) and we have taken ft, = 1. (Also, for 
the diffusion approximation to be valid, /cL 7 (t*) < 1.) 

The function I 2 (k) in the eqs.(|)-(|lg) is a dimension- 
less mode-coupling integral given by 



non-linear dependence of Ci on M{k) will imply a non- 
Gaussian statistics for the anisotropies. It should also be 
emphasised that in standard models the Cis have a sharp 
cut off for I > R^/Ls, due to Silk damping. But strong 
damping of Alfvcnic perturbations is expected only on 
scales smaller than VaLs q. Moreover, for small scale 
rotational perturbations, the damping due to the finite 
thickness of the LSS is also milder than for compressional 
modes. We note in passing that such rotational pertur- 
bations could also induce a "magnetic" type polarisation 
anisotropy. We will estimate this effect elsewhere. 

Several satellite and interfcrometric experiments are 
planned which will map the small angular scale 
anisotropies at the levels that these calculations predict. 
We have identified a new physical effect which can pro- 
duce a detectable signal in the microwave background 
on arc minute scales. If magnetic fields play a role in 
the development of large-scale structure in the Universe 
then we should be able to detect the effects of tangled 
primordial magnetic fields. 
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where |(k + q| = (k +q +2kq(i) ' . In general, I{k) can 
only be evaluated numerically. We will do this elsewhere, [2] 
but its order of magnitude contribution can be estimated 
from the following simple case. Suppose the magnetic 
spectrum has a single scale, with h(k) = fc<5o(fc — kg), 
where Sd{x) is the Dirac delta function). In that case 
< bg >— Bq and the field points in random directions 
but has a unique scale k^ 1 . The mode-coupling integral 
can be evaluated exactly for such a spectrum. We find 
thati(fc) = (fc//c )[l-(fc/2fco) 2 ] 1/2 , forfc < 2k , and zero 
for larger k. One can see that I(k) contributes a factor 
of order unity around k ~ k , with I(k ) = V3/2. For 
more complicated magnetic spectrum, with a multitude 
of scales, I(k) can be thought of as a superposition of 
these elementary contributions, and could be somewhat 
larger. 

From Eqs.(|s|)-(^o|), we see that in general for a tangled 
field of order B n - 3 x I0" 9 G, one expects a RMS CMBR 
anisotropy of order 5/iK or larger, depending on the con- 
tribution of I{k) and the value of /. The anisotropy in hot 
or cold spots could be several times larger, because the [5] 
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